CO-36-63-Soo»-a66(59)  MS 


CCMBIHATORIAL  METHODS  HJ  THE  THEOH?  OF  DAMS 

hr 

La Jos  Takace 

Columbia  University 
Department  of  Mathematical  Statistics 

March  25,  1963 


.  3  T  I  A 

RPR2  2 1963  I! 


y^j/Sli  w  Lb  I 

TIBIA 


This  research  was  supported  tgr  the  Office  of  Kaval  Research  under 
Contract  Humber  Horuv 266(59),  Project  Humber  0*2-205 „  Reproduction 
in  sbele  or  in  part  is  permitted  for  any  purpose  of  the  United 
States  OevenuMot, 


1.  IHTRGDOCTTOH 


In  this  paper  we  shall  be  concerned  with  tea  mathematical  models 
of  Infinite  dans.  In  the  first  model  independent  randan  inputs  occur 
at  regular  tine  intervals  and  in  the  second  model  independent  randoot 
Inputs  occur  In  accordance  tdth  a  Poisson  process  0  The  first  nodal 
has  already  been  studied  ty  J-  GartL  [61,  g.  F.  Teo  [161  and  others,  and 
the  second  model  by  J.  Gani  and  [73.  J»-Qe3.  aad  [8], 

D.  0.  Kendall  [93»  and  others.  For  both  models  ire  shall  find  explicit 
formulas  for  the  distribution  of  the  content  of  the  dan  and  that  of 
the  lengths  of  th©  wet  periods  and  day  periods.  The  proofs  are  ele¬ 
mentary  and  based  on  t&o  generalisations  of  the  classical  ballot  theorem. 


2.  OBKERAUZATIOSS  0?  TEB  CLASSICAL  BALLOT  THEOREM 


'Qio  classical  ballot  theorem  JLc  bo  fclltnsi  Suppose  that  in  a 
ballot  candidate  A  scores  a  vote'  sad  ooedidato  B  scores  b  votes. 
Let  a  >  pb  ifcero  \x  0  is  on  integer.  The  probability  that, 
throughout  the  oountingv  tho  number  of  votes  for  A  is  always  greater 
than  p  tinas  the  rasabar  of  votss  rogictsrad  fter  B  is 


P  B  &J*. 

r  a+b  • 


provide!  that  all  the  voting  records  ara  equally  probable, 

Fcrnala  (1)  far  p  »  1  was  found  in  1887  by  J.  Bertrand  [4]  and 
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tor  p  >  1  ,  also  la  1887,  by  £.  BSrMer  D3*  **»  proof  of  (1)  fbr 
p  »  1  w  given  In  1887  by  D,  Andrrf  [2]  and  for  p  >  1  in  1924  br 
A.  Aomill  [1% 

The  olaasioal  ballot  theorea  can  also  bo  fbr«latod  as  fellows i 

Lot  v  «  0  if  the  jvth  vote  la  east  for  A  and  let  v  *  (p+1) 
r  * 

If  the  iwth  vote  Is  cast  for  B.  1st  ns  a+b  and  k  »  b(vrf-l).  Then 

(2)  F  [  rr<x-tor  r  »  »!+...+  *n  »  *3  -  1  -  J  . 

if  0  <  k  <  n  . 

The  author  proved  by  nathsnaitleal  induction  that  (2)  also 
holds  if*  more  generally,  are  interchangeable  random 

variables  assuming  nonnogativo  Integer  values,,  ( 0'Z>  [10]  and  [U]«) 
Moreover  (2)  also  bolds  if  vn  are  cyclically  interchangeable 

random  variable b  assuming  nonnegative  integer  values.  For  this  latter 
oaee  a  simple  geometric  proof  -aas  given  ty  C.  _L„  yjallo^g  (oral  cam®!* 
nicotian).  Cf.  also  J„  C.  Tenner  [15]  and  ?V_D«&ag  [5]e  In  what 
feUaro  to  shall  prove  (2)  for  cyclically  interchangeable  random 
variables. 

theorem  i.  MjajgZBgBUagt  ‘i*  v*9  n  «&JBBBBags&a> 

lntogeiw valued  random  variables  and  that  all  the  n  cyclic  vemxtatlong 
of  (vI9  v2,  vn)  have  a  ooagpn  joint  dlStrltottOh.  Em 
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(3)  P  { vl+...«-vp  <  r  for  r  «  I,...,  nj  Vj+,.,+1^  »  k)  - 


1  -  £  il  0  <k  <n. 


rrrnM  flair  ffr  htrt  ilfr  to  frfW- 

PROOF,  Let  k^,  l^,,,,,  kQ  be  fixed  nonaegatire  integers  with 
son  k^*  kg*,..*  kn  ■  k  where  0<k<n.  Define  kj^  ■  kj  tor  J  *  1,2,..,  , 
life  shall  prove  that  among  the  n  cyclic  permutations  of  (kj,  V-  V 
there  are  exactly  ruk  for  which  the  sum  of  the  first  r  members  is 
less  than  r  for  all  r  =  1,2,, ,.,n„  Hence  the  theorem  lmaadiately 
follows  for  0<  k  <  n.  If  k  *  0  or  k  >n,  then  (3)  is  trivially  true. 

Define  A^  “  Xkj^  +  ...  +  k^)  for  J  *  1,2,,.,  ,  Then  ^3*n  “ 

A  ^  ♦  An  ■  Aj  +  Ca-k)  tor  j  “  1,2,..,  .  Let  m  be  the  greatest 
positive  integer  for  which  AR  ■  min(  A j,  A2,...,An).  Now  we  shall 
prove  that  there  are  exactly  twk  values  among  i  *  m+l,,..t  nrta  such 
that 

(4)  A  i  <  Aj  for  all  j  «  i+1,...,  i-m, 

that  is,  there  are  axaotly  n-k  permutations  among  (ki+1,...,k1+n)  , 

1  *  sH*l, o><)  m+n,  for  thich  the  sum  of  the  first  r  member*  is  less  than 
r  for  all  r  *  1,2,  •  9-9 

Denote  by  i^,  i2,...,i ^  ^e  greatest  indices  such  that 
-  Aa+i,  Aig  »  Aj  2,...,  «Aa  +  (n-k)  * 

respectively.  They  exist  because  -  A^  ■£  1  for  every  j.  Qr 

the  definition  of  a  we  have  1^  =  jrtn  and  therefore  m  <  i^  <  <,,,<  ln  >T 

hHo.  Clearly  (4)  holds  if  and  only  if  i  »  ij,  i2,..»,  ,  This 

proves  the  assertion, 

Now  we  shall  also  prove  a  further  generalisation  of  the  classical 


Allot  theorem*  The  following  theorem  wa  found  ty  the  author  [123 
fbr  interchangeable  random  variables,  however,  ve  shall  prove  it  here, 
'lightly  more  generally,  for  cyclically  interchangeable  random  variables. 
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the  symbol  [a]  denotes  the  greatest  integer  <  a.  Ry  Theorem  1 
we  hare  for  0  <  y  £  t  that 

(8)  1  -  $  £  P  }»{■>  +  ...  +v[‘)<rftrr«ll  ...f  2?  U<t)  «  yj  S 1 

because  1  «  1.2.  «...  2®,  are  cyclically  interchangeable  randon 

variables  that  assume  nonnegative  integer  values  and 

♦...+»« a  ^ 

if  X  (t)  »  y  o  If  m  — >  oo  in  (8),  thro  by  the  continuity  theorem 
of  probability  we  get  that 

(9)  P  {/*  (a)  |u  for  0  <  u  <  t  J/£(t)«yJ  *  1  - 

for  0  §y  «  This  proves  ($)  . 

REMARK,  in  [13]  we  proved  that  if  Vp  y^,  ....  are  inter¬ 
changeable  random  variables  that  assume  nonnegative  integer  values, 
then  the  probability  that  +  ...  +  <  r  holds  ftor  exactly  J  values 


among  r  *»  1,  ....  n  given  that  +  ...  +  vn  «  k  is 

(10)  Pj-£ 

P  i’l  *-*  Vl  -  1  1  *1  ♦•••♦  'n  '  “3 

iwuj 

if  0  <  k  <  n-1 

and  j  *  n-k,  ....  n-1,  and 

<u) 

p3  -  i 

if  k  «  n-1  and  j  *  1.2.  ....  n. 

Qy  using  the  same  proeedore  as  we  used  in  proving  Theorem  2  we  can 
obtain  from  (10)  and  (11)  the  following  result*  12  Z  v  Zv  ....  Zn 


■^m^i^fgSgs^p^r  > 
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are  nonnegative,  interchangeable  random  variables,  if  ^  (a)  is  defined 
hr  (6),  and  if  9  (t)  denotes  the  measure  of  the  set  [us  X (u)  <  u 
and  9  £  u  £  tj  ,  then 

(12)  P  (§>  (t)  £  x  l2T (t)  =  yj  •  ^  J  )  P  {u  </(u)  <  a-Hiu  J^Ct)  «  yj  dv 

t<U+7 

u£y,  v£x 

if  y  <  t  and  t-y  <  x  £  t,  and 

(13)  P  (^(t)  <x  |/(t)  -  t]  »f 

if  0  <  x  <  t  « 


3.  HEXHJLilR  INPUT 

Suppose  that  at  times  a  ■  1,  2,  ...  water  of  quantities  v...,  Vg,  ... 
is  flowing  into  a  dam  (reservoir)  and  the  release  is  continuous  at 
constant  unit  rate  when  the  dam  is  not  empty.  Suppose  that  r^,  Vg,  ...»  vQ,  . .. 
aro  identically  distributed,  mutually  independent  random  variables  that 
assume  nonnegative  integer  values.  Denote  ly  the  content  of  the  dam 
immediately  after  time  n,  Dion  we  have 

(14)  r(lnm  ^  7n-l  ”  *  vn  *  n  *  »•**  • 

The  initial  content  is  a  nonnegative  integer.  Denote  ty  0Q  the 


tine  of  the  first  emptiness,  i.e. ,  the  smallest  value  of  a  '-rich  tl'.'it 
r^n  »  0,  Allowing  the  initial  wet  period  (if  any)  dry  periods  and 
wet  periods  alternate.  Denote  by  Op  ©g,  ©n,  ■«  the  lengths 

of  the  successive  wet  periods  other  than  the  initial  one.  They  aro 
identically  distributed,  mutually  independent  random  variables.  The 
dzy  periods  are  also  identically  distributed,  mtually  independent  rardora 
variables  and  independent  of  the  wet  periods.  The  probability  that  » 
dry  period  has  length  k  is  [1-P  ~  Oj  ]  [P  [  »  0  3  3^”“  £°'~ 

k  1|  2)  wo  ®  # 

By  (14)  w©  have 

(15)  »  max  {vy  +  ...  +  vn-  (n-r)  for  r  =  1,  n  and 

"to  +  vl+  •••  +  vn~ni  ' 

In  what  follows  we  shall  use  the  notation  JSQ  »  0  and  Nn  “  \  +  »•«  *  \ 
for  n  *  1,  2,  ...  * 

THEOREM  3.  We  have 

■.  (  \  n**l  S^i”i  Q  r  1  ( 

(16)  pfyn  £  k  l^o  »  ij  -  p  *  n+k-ij  -  2_,  {1"  ^3)P 

j»l  -t*0 

and.  In  particular, 

f  {%>  - 0  ito  =  l3  *  £ (1  -  ii> p  &  -  j)  • 


(17) 


PROOF,  If  we  replace  y^  v2,  ...»  Yfl  by  vn,  v^,  ....  ^ 
respectively  la  (15),  then  we  obtain  a  new  random  variable 

(18)  *  aax  |llr  -  r*l  for  r  *  1,  n  and  Hn  -  n  +'^q j 
which  has  the  same  distribution  as  ^  T  •  Thus 

PhnSk  hoBi]  »  a+k-xj  - 

(19) 

P  [Kn<  n+k-i  and  Ny>  i**k  for  mim  r  *  1#  nj  „ 

Let  r  «  ;j  (J  *  1,  a«i)  be  the  greatest  r  for  i&ich  |  r+k<, 
Then  Nj  »  JHc  and 


(20) 


*iln$  k  l"(o  *  lj  ’  p  K  5  Wfc-ij  -  r  (»3  -  **j 

j»l 

P  <  n  -  j«.i  find  -  Nj  <  r  »  J  for  r  *  j+1 


By  Theoreia  1 

(21)  P  ;  B  «  H.  <  r  -  j  for  r  »  J+:.,  n  and 

s  *  J 

Hn  .  n3 s  w*ij  p  [V,  *(}  • 

£»0 


Patting  (21)  into  (20)  we  get  (16),  If,  in  particular,  k  *  0, 


than  'ay  Theorem  1 


(22)  P{^a*0|^0*ij  <  r*&  and  Ky  <  r  to r  r  *  X, 

£  P(Sn^  {1-5>  * 

jsaO 

THEOHEK  'l-o  SLteS^  2*£i 


23)  p  £  51  I  *  ij  »  2Z  d  P  :5  J,i) 

PECO?,.  Hoy  "bo  have 


P  [©^  <  n  j  a^o  a  *  ?  ^3^  *  zv?.  tor  soi so  ^  **  *j  «r<  o  ?*  ®  | 


?  fa-  *  5-l\  P  feu  S_  <  J-r  ftM?  i-  *  1,  ,  ,  „  5-1  j  SL  a  i«.:  ; 

k  ■"  - 

’Efi  Vcr  Th-soron  1  the  second  factor  in  the  .*.  sa  5.3  eqm?.  to  i/5* 


™»5.  HJ^Jar  5-  -  1.  2  ,,.  .  t&r; 


05)  *  {%  Snj  -p  -jli  ?  {»ja  »  3  I  %  »  °i  • 


PB01F,  It  can.  easily*  b«  aa-ra  iix:l 


P)<3^  <tx)  a  P{3^  <  3S-1  ior  se.e  ::  *  1;  0 


p  cVi  *  ;l  J  •%  *  °j  p  i3j+i  *  sv*>:  *  J'r  :*  “ 


and  (25)  fallow  free  Thaorwa  1* 


REMARK,  If  vo  suppose  -Quit  the  level  of  tho  dam  say  vary  in  the 
interval  (-  oo,  oo),  that  is,  tho  dam  novcr  beconas  empty,  than  the 
probability  that  la  the  tin*  interval  (0,  n)  the  total  tine  daring 
t*Jeh  tho  level  ie  below  the  initial  W  given  that  \  *  is 

(2?)  P  jHy  <  r  ibr  j  indices  r  «  1,  ,,,,  n  |  ^  «  kj  * 

> .  -a+SfiSiy ?  (Bi+i  ” 1  ^  \ 

If  0  <  k  <  n-1  and  j  «*  n-k,  nJ.,  and 

(28)  P^Hp<r  fbr  j  indices  r  *  1,  ...»  a  |  3^.  *  ».1$ 

If  k  *  n»l  and  J  B  I,  n*  These  follow  free  (10)  and  (11) 
respectively,, 


4.  POISSON  INPOT 

Suppose  that  in  the  tlae  interval  (0,  oo)  eater  is  flowing  into 
a  dan  (reservoir)  according  to  a  random  procase  end  the  release  is 
eontinaeos  at  constant  unit  rate  ehea  the  dan  ie  not  empty.  Denote 
tgr  Z  (t)  the  total  quantity  of  eater  flowing  into  the  dan  daring  the 
tine  interval  (0,  t).  It  is  supposed  that 


******  0  <  t  <  00}  is  a  Poisson  process  of  density  ^  and 

^ 2*  ****  •••  ***  identically  distributed,  actually  independent, 

positive  random  variables  and  independent  of  £v(t)  j  . 

Denote  by  ^(t)  the  contant  of  the  dam  at  ties  t0  The  initial 
content  ^(0)  =:  0.  Denote  ly  the  time  of  the  first  o?nptin»B8j  ©Q  =  0 
if  <^(0)  *  0„  SbllcfTsing  ths  Initial  wot  psvdod  (if  any)  dry  periods  and 
wat  periods  alternato.  Denote  by  6^,  ©2?  ®n,  ...  tha  lengths  of 

tha  saaoesalve  wet  periods  other  -than  the  initial  one.  Tfcay  are  identically 
distributed*  satnally  independent  random  variables.  Tho  dry  periods 
are  also  identically  distributed,  autuaTIy  indopendent  rendoea  varinblen 
and  independent  of  the  vet  periods.  Oho  probability  that  a  dry  period 
has  length  |x  is  X  «  e”^2  for  x>0„ 

First  we  nantion  that  hy  Theorem  Z  wa  have 

(30)  P  [X  (u)  £  u  Ibr  0$ult  IX(t)*yj  *  (1  - 

if  0|j  |t,  Par  if  we  know  that  in  the  time  interval  (0,t)  there 
are  n  >  0  events  in  the  Poisson  process,  then  the  oosurrenee  times 
^l*  ^2*  *9°*  ^n  hav*  ****  8aaQ  joint  distribution  as  the  coordinates 
arranged  in  increasing  order  of  n  points  distributed  uniformly  and 
Independently  of  each  other  in  the  interval  (0,t).  Thus  by  (5)  the 
probability  that  ^  (u)  £  u  far  0  £  t*  £  t  given  that  X  (t)  *  y  and 
v(t)  *  n  Is  (1-^)  if  0  §  y  £  tc  Since  this  probability  is  Independent 
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« f  the  oooditlon  r(t)  ■  n,  (30)  fellowe  iwwdiately. 

using  (30)  or  the  procedure  which  we  used  in  proving  Thao  ram  2 
»  obtain  the  following  theorems  corresponding  to  Theorems  3#  **■»  and  5* 
In  what  follows  we  shall  use  the  notation  P  £z(u)  <  xj  ■ 

P  f  x  <  £(u)  $  x+dx]  regardless  of  whether  u  depends  on  x  or  not* 

THEOREM  6*  c  >  0  ^  x^O,  then 

(31)  P  (  ^(t)  5  x  1  n^(0)  »  0]  «  P  <  t-*OMJ  j  - 

K  ^  ^  p [ZM  S a4xi  <V  p fcXt-u)  St}  , 

n+w<t-c 

Ogi,Ocv 


In  Particular. 

t-c 

(32)  P  h<t)  *  0  1  ^(0)  =  e j  *  ^  (1  -  ^  p  £*<t)  <  yj 

0 

If  t  >_  c  and  0  otherwise. 


THEOREM  7.  &  e  >  0,  then 

t 

(33)  p {«0  <  t  |^(0)  *  ej  -  ^  *  d^  P^(u)  <  iwoj 

« 

if  t  >  »ai  0  t  <  Co 
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THSOBZH  8.  Jl  t  >  0,  then 

t 

<3W  r  KstJ  •  i\  l4ttp{o<^(»)  juj 

0 

for  1  *  1,  Z9  ...t  and  X  =  lim  PfZ(t)  >  oj  /t. 

t->  0  L 

The  above  three  theorems  have  been  proved  in  [4]  under  the  more 
gensral  assumption  that  \X  (t)f  0  <  t  <  coj  is  a  stooisastic  procese 
with  noaneg6tives  stationary,  indepandent  increments, 

REMARK.  If  we  suppose  that  tho  level  of  the  den  may  vaiy  in  the 
interval  («oo,  oo ),  that  is,  the  dam  never  becomes  empty,  then  the 
probability  that  in  the  time  interval  (0,  t)  the  total  tine  during 
which  the  level  is  below  the  initial  level  given  that  X  (t)  «  y  is 
equal  to  (12)  or  (13)  where  now  [^(t),  0  <  t<  oo]  is  defined 
by  (29). 
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